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Abstract. The domain growth processes originating from noise-induced nonequilibrium phase transitions 
are analyzed, both for non-conserved and conserved dynamics. The existence of a dynamical scaling regime 
is established in the two cases, and the corresponding growth laws are determined. The resulting universal 
dynamical scaling scenarios are those of Allen-Cahn and Lifshitz-Slyozov, respectively. Additionally, the 
effect of noise sources on the behaviour of the pair correlation function at short distances is studied. 

PACS. 05.40.-a Fluctuation phenomena, random processes, noise, and Brownian motion - 64.60.-i General 
studies of phase transitions 



1 Introduction 

Growth processes for systems in their evolution towards a 
final state of thermodynamic equilibrium have been inten- 
sively studied during the last decades Q. A particularly 
interesting situation concerns the case in which equilib- 
rium corresponds to either of two equivalent phases of a 
scalar field 0(x, t). Since in equilibrium the relevant free 
energy F is minimized, "equivalent" stands in this context 
for phases that have the same (minimum) value of that 
free energy. A typical case corresponds to the evolution 
following a sudden quench from a very high-temperature, 
homogeneous phase, to a final state below a critical tem- 
perature for an order-disorder transition. In this process, 
one can distinguish three regimes: 

i) Right after the quench, the system finds itself out 
of equilibrium and small domains start to appear, corre- 
sponding to islands of each one of the two phases. The 
initial instability mechanism varies with the system. 

ii) At late times, motion of the domain walls separating 
equivalent phases leads to a situation of domain growth, 
also known as coarsening. This regime occurs when the 
typical domain size R{t) is much larger than the width of 
the interface between domains and much smaller than the 
system size. The main general mechanism for coarsening 
is the reduction of the excess of interfacial energy, which is 
accomplished by reducing the interface curvature. The ex- 
act details of the domain growth dynamics depend on the 
spatial dimension d and on the existence (or absence) of 
conservation laws. For dimension d > 1 and for a system 
whose order parameter is not conserved (i.e., which evolves 
towards a final one-phase, symmetry-breaking, equilib- 
rium ordered state), the domain size R(t) grows with time 
t as R(t) ~ t 1 !" 1 (Allen-Cahn law). On the other hand, 



for a system with a conserved order parameter (i.e., in 
which the final state has two coexisting phases), it has 
been found that R(t) ~ t 1 ^ 3 (Lifshitz-Slyozov law). More- 
over, in both cases the relaxation of these systems towards 
equilibrium has been found to be self-similar. This claim 
basically states that, after the initial transient regime, the 
only relevant length scale is the domain size R(t), so that 
two images of the system taken at two different times, t\ 
and £2, can be made to match (in a statistical sense) by 
rescaling their lengths by R(t\) and Rfa), respectively. 
The validity of this dynamical scaling law [g,[3) has a pre- 
cise statement in terms of the correlation function, and has 
been extensively verified in many experiments, as well as 
in computer simulations of discrete and continuous models 

iii) Finally, when the domain size R(t) is of the order 
of the system size, the system asymptotically reaches an 
equilibrium steady state, in which either only one of the 
two phases occupies the whole system (non-conserved or- 
der parameter case) or two large domains separated by a 
single boundary coexist (conserved order parameter case). 
In this case, the only relevant length scale of the problem 
is the equilibrium static correlation length. 

Quite generally, the growth processes mentioned above 
have been described theoretically by the following dynam- 
ical equation: 



dt 



= -r- 



5T 



(x,i) 



(1) 



where T[4>] is a suitable nonequilibrium free-energy func- 
tional, r is a kinetic coefficient and ?j(x, i) is a Gaussian 
white noise representing internal fluctuations, taken to 
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have zero mean and correlation: 

(r?(x, t) V (x', t')) = 2e <5(x - x') S(t - t>) . (2) 

When e = rksT (what is known as the fluctuation- 
dissipation relation), the system reaches an equilibrium 
state governed by the well-known Boltzmann-Gibbs dis- 
tribution exp[— J-/kB T]. In the absence of fluctuations, 
the free-energy T is a Lyapunov potential for the dynam- 
ical evolution, and the resulting dynamics is classified as 
potential Q . It is worth stressing that the growth laws and 
dynamical scaling introduced before hold independently 
of the strength of the fluctuation terms, i.e. of the value 
of e. This is true as far as those fluctuation terms are 
not extremely large (such that they inhibit the formation 
of order), or near a critical point (where fluctuations are 
largely amplified). 

Most of the work in the past has dealt with the sim- 
pler case when the corresponding deterministic dynamics 
is potential. Only in last years, attention has been focused 
on more general dynamics where the final state is not of 
thermodynamic equilibrium. Most studies have been con- 
cerned with the case in which the deterministic dynam- 
ics does not imply the minimization of a free energy, the 
so-called non-potential dynamics New and interesting 
effects can appear in this case. For instance, it has been 
found that several stable phases can coexist in the system 
in a non-coarsening situation Q. 

A new type of non-potential systems presenting order- 
disorder phase transitions has been introduced recently. 
The main feature of those systems is that the initial insta- 
bility leading to domain growth and coarsening is induced 
by fluctuation terms of external origin. This is a rather 
counterintuitive effect of noise, in that an increase of the 
noise intensity can lead to a spatially ordered pattern. This 
phenomenon is nowadays well established and has lead to 
a host of new effects in which noise has a spatial ordering 
role (see Q for a recent review). These effects include en- 
largement of the coexistence region in standard models of 
hase transitions |8| , pure noise-induced phase transitions 
noise-driven structures in pattern-formation processes 
, and noise-sustained pulse propagation JL1|, among 
otners. In all these cases, a nonequilibrium oraered pat- 
tern can emerge as a consequence of the fluctuating terms, 
and disappears when the noise source is switched off. The 
final state cannot be made to correspond to any known 
free-energy, and the equivalence of the different phases 
can be established, whenever possible, by symmetry argu- 
ments. These systems can be described by the following 
Langevin equation 



d<Kx,t) 
dt 



-r 



8T 



(x,t) 



+ *7(x,t), (3) 



where £(x, t) represents external fluctuation sources (mul- 
tiplicative noise, if g depends on 0), to be characterized 
statistically (see later). 

So far, studies of noise-induced spatial order have been 
concerned with the characterization of the final, nonequi- 
librium steady state (which would correspond to regime 



(iii) above). In particular, much effort has been devoted 
to determine the phase diagram and the critical properties 
of these systems, including the possible existence of new 
universality classes for the steady state. A great deal of 
our knowledge of the detailed behaviour of these systems 
comes from numerical simulations of the corresponding 
stochastic partial differential equations. On the theoreti- 
cal side, studies have used mean-field theories, dynamic 
renormalization group, and linear stability analyses [Q. 
Renormalization-group arguments, for instance, show that 
multiplicative noise is less relevant than additive noise, 
when both are present, to the universal behaviour, and 
thus it does not change the universality of the resulting 
critical points @. In this paper, on the other hand, we 
are concerned with the noise-induced dynamical evolution 
of both conserved and non-conserved systems towards a 
nonequilibrium steady state (which would correspond to 
the scaling regime (ii) defined above). We will be inter- 
ested in studying the growth laws and whether dynamical 
scaling holds in these cases. 

This paper is organized in the following way. In sec- 
tion II we study numerically the scaling properties and the 
growth of domains for the two-dimensional non- conserved 
Ginzburg-Landau model in presence of both external and 
internal fluctuations. In section III we characterize numer- 
ically the phase separation dynamics induced by external 
noise in the two-dimensional conserved Ginzburg-Landau 
model. Finally, section III is devoted to the study of the 
pair correlation function at short distances in the presence 
of noise sources. 



2 Non-conserved model 

As stated in the introduction, when a system in a high- 
temperature disordered equilibrium state is quenched be- 
low its order-disorder transition temperature, for which 
equivalent ordered phases can coexist, domains of the new 
equilibrium phases appear and grow. For those systems 
whose order parameter is not conserved, one of the do- 
mains grows until it fills the whole system, assumed fi- 
nite. The domain boundaries move with a translational 
velocity that has been found, for spatial dimension d > 2, 
to be basically proportional to the mean curvature of the 
boundaries, and independent of the free energy of the in- 
terface. This mechanism leads to the well-known Allen- 
Cahn law for the growth of the average linear size R(t) of 
phase domains: R(t) oc t 1 / 2 . The same law, with a differ- 
ent physical mechanism, has been shown to hold in other 
models as well [ [[p . 

A particular model for which the above results have 
been established is the non-conserved Ginzburg-Landau 
model (known as model A in the literature of critical 
phenomenajl^l). It is defined by Eq. (|l|) with a constant 
kinetic coefficient _T = 1, and with the following free en- 
ergy: 

1 

4< 



T= \ d: 



2^ 



(4) 
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For fixed values of the diffusion constant D and the in- 
ternal noise strength e, there is a critical value a c (D,e) 
such that for a < a c equilibrium corresponds to a sin- 
gle homogeneous phase (the disordered phase), whereas 
for a > a c two equivalent equilibrium phases (with the 
same minimum value of the free energy) coexist. Due to 
the symmetry of J-, the ordered phases are related by a 
global change of sign of the field <fi. The critical value a c 
vanishes for zero noise intensity a c (D, e = 0) = 0. 

A new ordering situation can arise if we let fluctuate 
the control parameter a — > a+£(x, t). Even when the mean 
value of a is still below the critical level (a) < a c , it might 
be possible, by choosing a sufficiently large intensity of the 
external fluctuations, to have a bistable stationary (but 
no longer equilibrium) distribution for the field. The two 
maxima of the stationary distribution are again symmetric 
under a global change of sign. 

In the presence of both internal and external fluctua- 
tions, the non-conserved Ginzburg-Landau model is ob- 
tained by substitution of the free energy functional (Q) in 
the Langevin equation (^). One finds 




Fig. 1. Spatial pattern for the stochastic model A with a 2 — 
0.4, e = 10" 4 , a = -0.2, and D — 1 at t=200. The square 
lattice has 256 x 256 cells of mesh size Ax = 1. The black and 
white domains correspond to symmetric phases. 



at 



acj) — cj) 3 + D V 2 4> + 0£(x, t) + 7j(x, t) , (5) Dij accounts for the discrete Laplacian operator 



where both additive and multiplicative noises are Gaus- 
sian, with zero mean. The additive-noise correlation is de- 
fined by Eq. (||). In general, the external multiplicative 
noise £(x, t) may have specific temporal and spatial cor- 
relation functions. For the moment, we will assume that: 



nn(i),j ^dSijj 



(8) 



(£(x, t) £(x', t')) = 2<j 2 c(x - x') 5{t - t') 



(6) 



where nn(i) represents the set of all the sites that are 
nearest neighbors of site i. The discrete noises rji(t) and 
£i(i) are still Gaussian with zero mean and correlations 



i.e. a delta-correlated function in time and a general cor- 
relation function c(x) in space. In this case, for given e 
and D, the critical value of a depends on a in such a way 
that its value is lowered with respect to the case with no 
external fluctuations, a = 0. In the limit of small e, a lin- 
ear stability analysis |l4|,[l5| and a mean-field type theory 
]l6|,p7| give the approximate result a c — — ct 2 c(0). Since 
ordered states appear for a > a c , by choosing parameter 
values satisfying — <r 2 c(0) < a < we have a situation in 
which nonequilibrium steady phases induced by external 
noise appear and grow (see Fig. |l|). If external fluctua- 
tions are switched off the ordered state disappears, since 
the (negative) value of a is now below the critical value. 
In this section, we want to elucidate if the growth of these 
domains still follows the Allen-Cahn law and if the system 
still exhibits a scaling regime. 

Most of the results of this section have been obtained 
by a numerical integration of the stochastic model defined 
by Eq. (^). In order to perform the numerical analysis, we 
redefine the model by considering a regular d-dimensional 
lattice with N d points and lattice spacing Ax: 



and 



dt 



a^-tf+DY, A, h + m (t) + & & (t) , (7) 



where fa (t) = <j)(xi , t) and only one index is used to label 
the cells, independently of the dimension of the lattice. 



( Vi {t)r lj {t'))=2ej^S{t-t') (9) 



&(t)t j (t'))=2a 2 j^L i 6(t-t>). (10) 

This last expression corresponds to an external noise with 
a white spatial correlation in the lattice, i.e. c(0) = 1 j Ax d . 
Since we will choose Ax = 1 in all the results shown in 
the paper, c(0) = 1 in what follows. 

After the discretization procedure described above, the 
resulting set of coupled stochastic differential equations for 
the variables fa, ... , 4>n has been integrated numerically. 
We have studied two cases of noise-induced non-conserved 
dynamics, corresponding to two different multiplicative 
noise intensities, a 2 — 0.4 and a 2 = 0.6. In both cases, 
the additive noise intensity is e = 10~ 4 , and a = —0.2. 
Since these values satisfy —a 2 < a < 0, the linear stability 
analysis indicates that the homogeneous state is unstable 
for these multiplicative noise intensities |l4|,|l5|]. We have 
compared these two cases with the purely deterministic 
model e = a = for a = 1, and with the stochastic model 
with only additive noise (er = 0) for a = 2 and e = 0.7. 
In all cases we have chosen the coupling coefficient to be 
D = 1. 

Our numerical analysis has been performed in a two- 
dimensional square lattice of 256 x 256 cells with periodic 
boundary conditions and mesh size Ax = 1, except for 
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the stochastic model with only additive noise, which has 
been considered in a lattice of 128 x 128 cells in order 
to compare with the results of Ref. Jll| . We have used a 
first-order Euler algorithm with a time step dt = 10 -2 
in the deterministic case, and dt = 5 • 10~ 3 in presence 
of noise sources. The noises have been generated using a 
numerical inversion method Jl9| . As initial condition, we 
have chosen in all cases a random uniform distribution 
of the field in the interval (—1, 1) in order to simulate a 
high-temperature one-phase state. Our results are aver- 
aged over 15 samples in the deterministic case, 40 samples 
in the purely additive-noise case, and 30 samples in the 
two multiplicative-noise cases. 

In order to study the growth process and to examine 
the existence of dynamical scaling, we define the corre- 
lation and structure functions. These two functions give 
information about the spatial structure of the system at a 
certain time. By knowing how these functions evolve, one 
can obtain information about the growth of domains, and 
check if there is a scaling regime. In our discrete space, we 
define the pair correlation function as: 



(11) 



where the brackets denote an average over different initial 
conditions and realizations of the noises. 

The structure function is the (discrete) Fourier trans- 
form of the pair correlation function: 



S(k fi ,t) = (Ax) d J2e- ir ^ G{r . jt) 



(12) 



where /i labels the N d points in Fourier space. 

In practice, the structure function is computed in an 
easier way using the equivalent definition 



:(|<KW)| 2 }, 



Ax d N d 

in terms of the Fourier transform of the field: 



(13) 



(14) 



and the correlation function G(ji,t) is then computed as 
the inverse Fourier transform of 5(k M , t). 

We perform a spherical average of the pair correlation 
and structure functions 



r<ri<r-\-Ar 



S(k,t) 



N, 



E s ow)' 



(15) 



(16) 



k<k^<k+Ak 



where the sums run over the set of lattice points (N r and 
Nk) between two circles of radius r and r + Ar, or k and 
k + Ak, in real and reciprocal space, respectively. 



Finally, in order to compare results from different pa- 
rameter values, we have used the following normalization 



9(r,t) 



G(r,t) 
G(0,t) 



s(k,t) 



S(k,t) 
G(0,t) 



(17) 



When the interface between domains is very thin com- 
pared with their size, the system has a single character- 
istic length R(t), which is related to the average size of 
domains. The scaling hypothesis states that the pair cor- 
relation function depends on time only through the time- 
dependent characteristic length: 



g(r,t)=g(r/R(t)) 



(18) 



The length R(t) can be defined in several ways, but in 
the scaling regime all of them should lead to the same law 
for domain growth. We have chosen R(t) as the distance 
at which the pair correlation function is half its maximum 
value, i.e. g(R(t), t) = 1/2. Since the structure factor is the 
Fourier transform of the correlation function, we derive its 
scaling law as 



s(k,t) = R(t) d s{kR(t)) 
with no other explicit dependence on time. 



(19) 
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Fig. 2. Pair correlation function for the 2-dimensional non- 
conserved Ginzburg-Landau model in the deterministic and 
stochastic cases. The dotted (t=300) and dashed (t=500) lines 
correspond to the deterministic case with a = 1 and D = 1. 
The solid (t=300) and the long-dashed (t=500) lines corre- 
spond to the stochastic model with e = 10~ 4 , a 2 = 0.6, 
a = —0.2 and D = 1. In presence of noise sources the pair 
correlation function exhibits a peak at r = 0, as shown more 
clearly in the inset. 

We plot in Fig. the pair correlation function for the 
deterministic case with a = 1 and for the nonequilibrium 
stochastic case with a 2 — 0.6, a = —0.2, for two differ- 
ent times. As time increases, the correlation fuction tends 
to a monotonically decreasing function with an increasing 
characteristic decay length. A closer look (see inset of Fig. 
||) shows that in the stochastic case the pair correlation 
function is not smooth near the origin. This behaviour 
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Fig. 3. The Allen-Cahn law. Stars correspond to the deter- 
ministic model with a = 1 and d = 1. Crosses correspond to the 
additive case (e = 0.7, a — 2 and D = 1), squares to a 2 — 0.4, 
and triangles to a 2 = 0.6. The two cases with multiplicative 
noise have e = 10~ 4 , a = —0.2 and D = 1. The solid line is a 
guide to the eye. 

at short distances appears whenever noise sources (either 
multiplicative or additive) exist, and it differs from what 
is observed in the deterministic case, where the pair cor- 
relation function is smooth at all distances. As we will 
see, this behaviour is a bulk feature and, since it does not 
present scaling, we have avoided it in the following sec- 
tion by substituting G(r,t) near r = by a parabolic fit. 
The last section is devoted to the understanding of the be- 
haviour of the pair correlation function at short distances 
in the presence of noise sources. 

Figure |^ shows the square root dependence on time of 
the average domain size. The Allen-Cahn law is seen to be 
verified by all models under study. The prcfactor seems to 
depend on the noise intensity. 

Numerical results for the scaled pair correlation and 
the structure functions are shown in Figs. [|(a) and|I](b), 
respectively. Our numerical simulations verify scaling for 
both the pair correlation and the structure functions. Fur- 
thermore, the observed scaled structure functions agree 
with the theoretical prediction of Ref . [^o) , which has two 
adjustable parameters (a scale factor in each axis), and 
are also in agreement with Ref. [Il8| . 
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Fig. 4. Scaled pair correlation function (a) and structure fun- 
tion (b) for the non-conserved Ginzburg-Landau model. Bro- 
ken lines correspond to the deterministic model with a = 1 
and D = 1 for t=1000 (dashed line) and t=1500 (dot-dashed 
line). Empty symbols, squares (t=500) and circles (t=600), 
correspond to the additive model (e = 0.7, a = 2 and D = 1). 
Full symbols correspond to the nonequilibrium models. Squares 
(t=800) and circles (t=1000) correspond to a 2 = 0.4, while tri- 
angles (t=600) and diamonds (t=1000) correspond to a 2 = 0.6. 
Both cases have a — —0.2, D = 1 and e = 10 -4 . The continu- 
ous line in (b) is the theoretical prediction of Ohta et al. 



3 Conserved Model 

The conserved version of the model studied in the previous 
section, Eq. (0), is 

= -V 2 [o0 - 3 + D V 2 + <P e(x, t)] + „(x, t) , 

(20) 

which in the absence of external noise is the Cahn-Hilliard- 
Cook model pl[ and in the literature of critical phenom- 
ena |13| is known as model B. The conservation law is 
V^ 1 JJ7<ix0(x, t) — 4>o, constant, where V is the total vol- 
ume of the system. This constant is given by the initial 



volume fraction of the system. This conserved model is 
suitable for the description of the evolution of systems 
such as binary alloys, in which the total concentration of 
each component of the alloy is kept constant. In this case, 
the local order parameter, </>(x, t) represents the local dif- 
ference of concentrations of each component of the alloy. 
As before, we consider external fluctuations on the control 
parameter a. Both additive and multiplicative noises are 
Gaussian, with zero mean and correlation given by 

(r/(x, t) n(x, 0) = -2 £ V 2 S(x - x') 6{t - t') . (21) 
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and Eq. (|6|) . The Laplacian term in Eq. ( |21| ) ensures that 
in absence of multiplicative noise and if e = fc^T (the 
fluctuation-dissipation theorem), the stationary distribu- 
tion is still governed by a (restricted) Boltzmann-Gibbs 
distribution exp[— T/Ub T](5(y~ _1 J v dx(p(x, t) — <j> ). Let 
us consider this equilibrium conserved model in a high- 
temperature, disordered, one-phase state corresponding 
to a sufficiently low value of the parameter a. If we now 
suddenly decrease the temperature of the system below its 
transition value (equivalently, by increasing a), the disor- 
dered one-phase state becomes unstable and the system 
develops domains of the new phases which slowly tend to 
the equilibrium ones. The classical theory of the kinetics of 
first order phase transitions distinguishes two regimes for 
the initial instability leading to coarsening: for a quench 
taking the system deep inside the coexistence curve (crit- 
ical quench or small order parameter fa) the system is 
unstable against long-wavelength perturbations and un- 
dergoes a process of spinodal decomposition. On the other 
hand, for quenches close to the coexistence curve, the sys- 
tem evolves by nucleation and growth of the nuclei formed. 
In any case, since the order parameter is conserved, the 
growth of a given domain is at expense of the smaller do- 
mains of the same phase. Domain growth appears as a 
result of diffusion accross the interface between domains, 
caused by the interface curvature. The equilibrium final 
state is a two-phase coexisting state. In this time regime, 
the average domain size grows with time as R(t) cx i 1 / 3 , 
which is called the Lifshitz-Slyozov law. This power-law 
behaviour has been derived analytically for small volume 
fractions, although it has been seen numerically to be also 
satisfied even for large volume fractions 22 . As in the 



non-conserved case, when the growth of domains verifies 
the corresponding power-law behaviour, the pair correla- 
tion function and the structure function satisfy the scaling 
hypothesis defined in Eq. (|l^). The scaled functions have 
not been analytically obtained so far. 



In this section we are concerned with the dynamics 
towards the nonequilibrium, two-phase coexisting steady 
state induced by external fluctuations. Let us now con- 
sider the above nonequilibrium model [Eq. (p0[)] in a dis- 
ordered, 4>q = 0, one-phase steady state corresponding to 
a small intensity of the multiplicative noise. If we now in- 
crease abruptly the intensity of the external fluctuations, 
a, the system develops domains of new phases that grow 
and tend to the new stationary values (see Fig. ||). For 
an external noise with a white spatial correlation in the 
lattice (c(0) = 1, as in the previous section), a linear sta- 
bility analysis |23 shows that order appears in this case 
for values for a above a critical point a c = —2d a 2 . Hence, 
for —2d a 2 < a < 0, phase separation induced by fluctu- 
ations appears. We will now address the issue of whether 
this dynamics still verifies the Lifshitz-Slyozov law, and 
whether a scaling behaviour for the pair correlation and 
structure functions exists. 




Fig. 5. Spatial pattern for the stochastic model B at t=10000, 
with multiplicative noise intensity a 2 = 0.1, e = 10~ 4 , a = 
—0.2, D = 1 and <f)o — 0. The square lattice has 256 x 256 cells 
of mesh size Ax — 1. The black and white domains correspond 
to symmetric phases. 



As in the previous section, we define the model [given 
by Eq. (p0|)] in a discrete space of mesh size Ax 



dt 



k 



D 



ik 



a4>k - <j>% +D >^ Dkj 4>j- 



4>k £,k(t) 



rik{t) (22) 



The discrete noises r\i (t) and £j (t) are still Gaussian, with 
zero mean and correlations given by 



(Vi(t)Vj(t')) 







(23) 



and Eq. (|Io|). These conserved noise term 7?j(t) can be 
generated as the divergence of a stochastic vector field 

We have studied two cases of noise-induced conserved 
dynamics corresponding to two different intensities of the 
multiplicative noise, a 2 = 0.1 and a 2 — 0.2. In both cases, 
additive noise is e = 10 -4 and the control parameter is 
a = —0.2, which corresponds to a situation in which the 
disordeded state, 4>i — Vi, is the deterministically sta- 
ble phase. However, since —4a 2 < a < 0, it turns out 
that for these multiplicative noise intensities, the one- 



phase disordered state becomes unstable 15 . We have 
compared their behaviour with three equilibrium cases: 
the purely deterministic model for a — 2 and a = 0.2, and 
the stochastic model B with only additive noise (a = 0, 
a = 2 and e = 0.7). In all cases we have taken the coupling 
coefficient as D = 1. 

We have performed numerical simulations of the full 
model in a two dimensional regular square lattice of 256 x 
256 points with periodic boundary conditions and mesh 
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Fig. 6. Normalised pair correlation function for the 2- 
dimensional conserved model in the deterministic and nonequi- 
librium stochastic cases. The dotted (t=7000) and dashed 
(t=10000) lines correspond to the deterministic case with a — 2 
and D = 1. The solid (t=7000) and the long-dashed (t=10000) 
lines correspond to the stochastic model with e = 10~ 4 , 
a 2 = 0.1, a = —0.2 and D = 1. In presence of noise sources, the 
pair correlation function is not smooth at short distances, as 
seen more clearly in the inset. In both cases we have considered 
0o = 0. 

size Ax — 1. Since the model is self-averaging and the 
size we have chosen is relatively large, it has not been 
necessary to consider many different realizations in or- 
der to have small statistical errors. In particular, we have 
calculated averages over up to 6 realizations for the equi- 
librium models, and over 10 realizations for the nonequi- 
librium stochastic models with multiplicative noise. We 
have implemented a first-order Euler algorithm and used 
a time step of dt = 10 -2 in the deterministic case and 
dt = 5 • 10 -3 in presence of noise sources. As initial con- 
ditions we have considered a Gaussian random field with 
zero mean and variance 10 -2 (<fo = 0). 

The definition of the pair correlation function and the 
structure factor are the same as in the non-conserved case. 
The only difference is that now the spatial average of 
</>(x, t) is a constant determined by the initial condition. 
The pair-correlation function oscillates and decays to zero 
with a characteristic distance that is a measure of the typ- 
ical domain size. This is shown in Fig. FJ, where we plot 
the evolution of the pair correlation function for the de- 
terministic (a = 0.2) and stochastic (e = 10~ 4 , a 2 = 0.1, 
a = —0.2) cases. The first zero crossing defines the growth 
length, g(R(t),t) — 0. As in the non-conserved case, we 
have also found the same behaviour of the correlation 
function at short distances that differs from what is ob- 
served in the deterministic case (see inset of Fig. FJ). We 
have proceeded in a similar way as before by using a 
parabolic fit of G(r, t) at r = 0. 

As shown in Fig. ^, our numerical results indicate that 
the average domain size obeys the equilibrium Lifshitz- 
Slyozov law even in the presence of multiplicative noise. 
The corresponding scaled pair correlation and structure 



Fig. 7. Average domain size R(t) for the conserved model. 
Stars correspond to the deterministic case with a = 2 and 
empty circles to the stochastic equilibrium model with a — 2 
and e = 0.7. Full squares correspond to the nonequilibrium 
stochastic model with a 2 = 0.2, e = 10~ 4 and a — —0.2. The 
solid line is a guide to the eye. In all cases the Lifshitz-Slyozov 
power law behaviour seems to be satisfied. 

functions are plotted in Figs. ||(a) and ||(b). The scal- 
ing behaviour is seen to be verified by the nonequilibrium 
model with a 2 = 0.2. Furthermore, the scaled functions 
agree with the equilibrium ones. 

4 Noise Effects on the Short- Distance 
Behaviour of the Correlation Function 

As we have seen in the previous sections, the discrete two- 
point correlation function of systems undergoing domain 
growth in the presence of fluctuations is not smooth at 
short distances, for both conserved and non-conserved dy- 
namics (Figs. U and|^). This generic feature can be found 
in the existing literature [ p5| , but as far as we know, it has 
not been explained yet. As we will see in this section, this 
behaviour is mainly a bulk feature, and the interfaces do 
not play any role. Therefore, it can be studied considering 
the bulk stationary situation. 

4.1 Non-conserved Model 

We intend to show that the presence of fluctuations in 
the steady bulk solution is responsible for the non-smooth 
behaviour of the discrete correlation function. To that end, 
we linearize the model given by Eq. (FJ) around the bulk 
phase, which, in order to simplify the analysis, we take as 
the homogeneous null phase (0 = 0). Let us start with 
the non-conserved model with only additive noise, and 
linearize it around = 

^g-^ a 0+/JV 2 + r?(x,t), (24) 

where the additive noise has zero mean and correlations 
given by Eq. (FJ) . The stability of this model requires a to 



8 



M. Ibanes et al.: Dynamics and Scaling of Noise-Induced Domain Growth 




(b) 



Fig. 8. Scaled pair correlation function (a) and scaled struc- 
ture function (b) for the conserved model. Lines correspond to 
the deterministic model with a = 2 for t — 1000 (continuous 
line) and t — 1600 (dashed line). Empty symbols correspond to 
the equilibrium stochastic model with additive noise (a — 2 and 
e = 0.7) for t = 2000 (circles) and t = 3000 (squares). Full sym- 
bols correspond to the nonequilibrium stochastic model with 
a 2 = 0.2, e = 10" 4 and a = -0.2 for t = 7000 (triangles) and 
t = 9000 (diamonds). In all cases D = 1 and 4>o = 0. 



be negative. We now look for the dynamical equation of 
the pair-correlation function G(r, i), defined as the con- 
tinuum version of ( |TT|) (r = x — x') |Q: 



dG{r,t) 



2aG(r,i) + 2£>V 2 G(r,i) + 2e<5(r), (25) 



which has a delta contribution due to additive noise. In the 
stationary state, this equation can be solved in a straight- 
forward way in Fourier space. The resulting integral is 
divergent in d = 2, indicating the absence of a contin- 
uum limit in this case. On the other hand, the integral 
converges in d — 1. Numerical simulations of Eq. ( |24| ) [see 
Fig. ^(a)] show a strong dependence on the noise intensity 
of the slope of the correlation function at r = 0. Therefore, 
this is a non-scaling feature, which can be characterized by 
studying the spatial derivative of G(r, t) near the origin. 



From (Eq), we obtain 



dr 



G(r,t) 



r=0+ 



e 

2D ' 



(26) 



where we have taken into account that the pair correlation 
function is symmetric, G(r, t) = G(— r, t). This expression 
shows clearly the relevance of the noise intensity e. 

After this back-of-the-envelope calculation, we perform 
a more exhaustive analysis in discrete space, in order to 
compare it with the numerical results shown in Fig. ^|(a). 
We now consider both additive and multiplicative noise 
sources and take the bulk phase = 0. Thus, we linearize 
Eq. (0) around this phase and obtain 



dt 



(27) 



From this equation, and making use of Novikov's theorem 
, we can obtain the dynamical equation for the discrete 
correlation (4>i(t)(j)j(t)): 



d_ 

dt 



k 



Djkifrfa)) + 2a 2 <0i0j)- 



Ax d 



2e- Sij 



~Ax d 



(28) 



Since we are analyzing a bulk feature, the dynamics is no 
longer relevant, and we can perform our calculations in 
steady state. Therefore, we set ^{<j>i{t)<i>j{t)) = 0, for any 
i and j. Since we are interested on very short distances, 
we focus on i = j, to obtain: 



— y 

Ax 2 V x 

nn(i) 



((4>i^nn(i)) ~ (4>i)) + 



Ax d 



?> + 



Ax d 



= 0. 



(29) 



where the Laplacian term has been made explicit and the 
sum runs over the 2d nearest neighbors of i. 

In discrete space, the first derivative of the two-point 
correlation function Gi = G(r^) at i — is 



Gq 



(4>i^nn(i)) ~ (4>i) 



Ax 



(30) 



Taking into account the isotropy of the system, we have 
(<l>i<f>nn(i)) = G i i V nn(i) (31) 

where Gi is the two-point correlation function for dis- 
tances equal to one cell. 

Considering Eqs. (30) and (|3T|), expression ( p9| ) can be 
rewritten in the following form 



Gn — — - 



2dD 



-Ax Gn - 



dDAx 



Gn — ■ 



2dDAx 



d-l 



(32) 



If we now take the continuous limit Ax — > and consider 
that only additive noise is present (tr 2 = 0), we recover the 
continuous result (||^) given above in the case of d = 1 . 
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The numerical simulations shown in Fig. ^(a) have 
been performed on a one-dimensional lattice of N = 16384 
points and mesh size Ax — 1, for a = — 2 and D = 1. We 
have used a first-order Euler scheme with dt = 5 ■ 10~ 3 , 
and computed the two-point correlation function in the 
stationary state. 



G(r) 



0.0 1.0 2.0 3.0 4.0 5.0 

r 



(a) 



G(r) 



(b) 

Fig. 9. Stationary pair correlation function for the non- 
conserved (a) and the conserved (b) linear models for sev- 
eral noise intensities: only additive noise with e = 0.1 (solid 
line), e = 0.4 (dotted line) and e = 0.7 (dashed line), additive 
(e = 0.1) and multiplicative noises with a 2 = 0.2 (long dashed 
line) and a 2 — 0.5 (dotted-dashed line). In all cases a = —2 
and D — 1 and for the conserved case (j>o = 0. 



Figure hd shows numerical and theoretical results for 



4.2 Conserved model 

Following the same scheme as in the previous section, we 
now analyze the same phenomena for the conserved model. 
By linearizing Eq. ( |20|) with only additive noise around 
4>(r, t) = 0, we obtain 



d^(x,t) 
dt 



-V 2 [ac/)+ DV 2 4>] +r?(x,t), 



(33) 



where a must be negative and the additive noise is given 
by Eq. (|2l|). The corresponding dynamical equation for 
the pair correlation function is 



dG(r,t) 
dt~ 



-V 2 [2aG(r,t) + 2DV 2 G(r,t) + 

2eS(r)}. (34) 

In the stationary state we have 

aG{r,t)+ D\7 2 G{r,t) + eS{r) = h (35) 

where h is a constant since the steady state is homoge- 
neous. In the case cf>o = 0, this constant is zero M , B6[ . 
Comparing Eq. (3^) with Eq. ( |25| ) in the stationary case, 
it can be seen that also in the conserved case Eq. ( |2^ ) 
holds in d = 1 . 

In a more detailed analysis, we consider the lattice 
case when multiplicative noise is also present. We linearize 
Eq. (p2h around 4> = to obtain: 



dt 



^ Dik ( a <j) k +D ^2 D kj fa + & £ fc (t) 

+m{t), (36) 

where the discrete noises r\i (t) and £i (t) are Gaussian with 



zero mean and correlation given by Eqs.(|23j) and (10) 

The procedure of the analysis is the same as in the non- 
conserved case. First we look for the dynamical equation of 
the two-point correlation function, which using Novikov's 
theorem can be derived to be: 



-a{4>j(j) s )-D y ^ j D s u{(t>k4>j)- 



a2 ^2^d ( D sk(<ftk<f>j) + D jk {(f) k <p s ) 



Gq = (Gi—Gq)/Ax in two different cases, the linear model ^2 "sfc ( jj i± aaiq. ia m \ 
with onlv additive noise and with both a.dditive and mill- 2- -< Ax d ^ 8 ' S 



with only additive noise and with both additive and mul- 
tiplicative noises. Empty symbols are simulation results: 
squares correspond to different additive noise intensities 
(bottom axis) and a 2 = 0, and triangles correspond to dif- 
ferent multiplicative noise intensities (top axis) and fixed 
additive noise e = 0.1. Lines are the theoretical prediction 
(Eq. (|3^)). We have to note that our theoretical prediction 
uses the simulation result for Gq. 



Ax d 



(37) 



We now focus on the stationary state ■§i(4>i4 l j) = an d 
consider the case i = j 



2^A 



a((f>i(f) s ) + D} j D sk (4> k (t)i) +e 



Ax d 

-0, (38) 
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where the additive noise term has been introduced inside 
the Laplacian. This equation is satisfied if each term inside 
the Laplacian is equal to a constant h. 



- 2 E 



a{<j>i4> s ) + D^Dsk^kfa) + 

k 



Ax d 



DskiMi) + D ik (</> k <l> a ) )=h (39) 



for any s, i. As before, for cf>o — 0, the constant h must be 
zero. We now take the case s = i. If we use definition ( |30| ) 
for the discrete derivative of the pair correlation function 
and take into account the isotropy of the system, Eq. (pl[), 
we can obtain from Eq. (p9|) 



6*o 



aG Ax 2a 2 G Q e hAx 

2dD ~ DAx d + l ~ 2dDAx d - 1 ~ 2dD 



(40) 



which in the continuous limit Ax — > and in absence 
of external fluctuations (a 2 — 0) coincides with the con- 
tinuous analysis of Eq. ( p6| ) for d = 1. Notice that the 
dependence on the additive noise intensity is the same 
as in the non-conserved model, while the effect of multi- 
plicative noise is different. Both results can be made to 
coincide if the intensity of the multiplicative noise of the 
non-conserved model A, a 2 A is related to the one of the con- 
served model B, a 2 B , by the relation a 2 A = (2d/Ax 2 )a 2 3 . 
This relation between both models in presence of multi- 
plicative noise has been established in several analyses jl5], 

Again, we have performed simulations of Eq. (p6J) with 
a = —2 and D = 1 for <po = with only additive noise, and 
with both white additive and multiplicative noises in a one 
dimensional lattice of mesh size Ax = 1 and N — 16384 
cells. We have used a first-order Euler algorithm with time 
step dt = 5 ■ 1CP 3 . Results for different cases of the sta- 
tionary correlation function are shown in Fig. ^|(b). When 
only additive noise is present, the pair correlation function 
of both conserved and non-conserved models coincide [cf. 
Fig. H(a)]. As in the non-conserved case, the contribution 
of additive noise is here more important than the corre- 
sponding to multiplicative noise. 

Analytical and numerical results for G' are shown in 
Fig. [h]. Full symbols are the simulation results in this 
case: squares correspond again to different additive noise 
intensities (bottom axis) with a 2 = 0, and triangles to dif- 
ferent multiplicative noise intensities (top axis) with fixed 
additive noise e = 0.1. Lines are the corresponding theo- 
retical predictions [Eq. (^) with h = 0] using the values 
of Gq obtained from the simulations. 
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Fig. 10. Discrete first derivative G' r of the pair correlation 
function in the stationary state at r = for both non-conserved 
and conserved models with a — —2 and D = 1. Symbols are 
simulation results for both the non-conserved (empty symbols) 
and the conserved models (full symbols) and lines are the the- 
oretical predictions. Squares correspond to different additive 
noise intensities e (bottom axis) with a 2 = and triangles to 
different multiplicative noise intensities a 2 (top axis) with fixed 
additive noise e = 0.1. For the conserved model we have taken 
<t>o = 0. 



size R(t) of the domains, our results indicate that R(t) 
grows with time as a power law. The exponent of the power 
law is 1/2 for the non-conserved case and 1/3 for the con- 
served case. These two values coincide with the equivalent 
ones obtained in the decay towards an equilibrium state 
(either deterministic or in the presence of internal fluc- 
tuations). We have found a scaling description consistent 
with the fact that R(t) is the only relevant length scale of 
the system, similarly also to the equilibrium cases. More- 
over, the corresponding scaling functions are in agreement 
with the ones observed in the decay towards equilibrium. 
These results indicate that the physical mechanisms un- 
derlying noise-induced phase-ordering processes are iden- 
tical to the deterministic ones, namely interface curvature 
in the non-conserved case, with the addition of diffusion 
in the conserved situation p7| . Finally, we have analyzed 
in detail a non-smooth, non-scaling behaviour induced by 
fluctuations that appears in the pair correlation function 
at short distances. We have shown that this behaviour is 
due to fluctuations in the bulk, and provided a theoretical 
expression for its magnitude. 



5 Conclusions 

In this paper we have studied the dynamics of phase or- 
dering that appears in a non-potential situation induced 
by external noise. We have considered both non-conserved 
(order-disorder) and conserved (phase-separation) situa- 
tions. Concerning the time evolution of the characteristic 
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